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Motivated by recent experimental measurements, we study the quasiparticle spectra and the 
collective excitations in doped Si^IrCH, in which the interesting interplay between the electronic 
correlations and strong spin-orbital coupling (SOC) exists. To include the SOC, we use the Hugen- 
holtz diagrams to extend the fluctuation exchange (FLEX) approach to the case where the SU(2) 
symmetry can be broken. By using this generalized FLEX method, we find a weak pseudogap be¬ 
havior near ( 7 r, 0) in the slightly electron-doped system, with the corresponding Fermi arc formed 
by the partial destruction of Fermi surface. Similar features also appear in the hole-doped system, 
however, the position of the Fermi arc is rotated 45° with respect to the former. These results are 
consistent with the recent angle-resolved photoemission spectra in Sr 2 lr 04 . We elaborate that these 
anomalous phenomena are caused by the scatterings of quasiparticles off the isospin fluctuation 
derived from the effective J e fr = 1/2 doublet. 

PACS numbers: 71.27.+a, 71.70.Ej, 71.10.-w, 71.18.+y 


I. INTRODUCTION 

Recently, the 5 d transition-metal iridium oxides have 
attracted significant attention, because they exhibit a 
number of exotic phenomena induced by the spin-orbital 
coupling (SOC) and the correlation effects of electrons^. 
Of these iridates, Sr2lr04 of particular interest for 
it shares many analogies with the parent compound 
La2Cu04 of high- T c cuprate superconductors, such as the 
same layered perovskite crystal structure^, the same anti¬ 
ferromagnetic (AFM) insulating ground stated, the sim¬ 
ilar magnetic excitation spectrum^ and electronic struc¬ 
ture^. It has also been proposed theoretically to realize 
the unconventional superconductivity via doping^. 

In particular, the recent angle-resolved-photoemission- 
spectroscopy (ARPES) measurements on the slightly 
electron-doped S^IrCU show a temperature-dependent 
pseudogap phenomenon^. The intensity of the spectra 
is much suppressed in an extended region near (0,7r), 
resulting in the Fermi arc which resembles the case in 
the hole-doped cuprates^. For the slightly hole-doped 
S^IrCU, another ARPES experiment also exhibits the 
trait of pseudogap^, but with the suppression of spec¬ 
tral near the (7t/2,7t/2) point. As the pseudogap puzzle 
is the longstanding unsolved problem in cuprate super¬ 
conductors, studying the origin of the Fermi arcs and 
pseudogaps in doped Sr2lr04 is not only important for 
this material itself, but might also help to investigate 
the similar phenomena in high-X/ cuprates and/or other 
materials. 

Sr2lr04 shows a multi-orbital electronic structure 
where the t 2 g and e g orbitals are separeted by large crys¬ 
tal field. The five electrons (one hole) reside in the lower 
t 2 g manifold of xy, xz , yz orbitals. In spite of the large 
band width and small Coulomb interactions, S^IrCU is 
an antiferromagnetic insulator— A It has been proposed 
that the SOC breaks this sixfold degenerate manifold 


into completely filled J e g = 3/2 bands and a half-filled 
Jeff =1/2 band (Kramers doublet) which is further split 
by the relatively small Coulomb interactions. Thus, it 
can be simplified to an effective one-band half-filled sys¬ 
tem, which hosts an isospin J e g =1/2 spin-orbital Mott 
insulating ground stat e — . However, the validity of this 
isospin J e ff =1/2 Mott picture including its robustness 
to dopings still remains an open question^—. 

Motivated by these experimental and theoretical 
progress, we carry out a theoretical study of the collec¬ 
tive excitations and quasiparticle spectra in the doped 
S^IrCU, based on the three-orbital Hubbard model with 
the inclusion of the SOC. 

The multi-orbital structure together with the SOC 
confines the exact diagonalization and quantum Monte 
Carlo (QMC) methods to very small systems. In view of 
this, the fluctuation-exchange (FLEX) approximation is 
a good alternative^—. The FLEX method has advan¬ 
tages to handle various collective fluctuations, and the 
calculation results agree well with the QMC simulations 
for the Hubbard model with the moderate on-site inter¬ 
action I / 17 ’ 18 . Up to now, the FLEX approach has been 
extensively applied to the high-T/ cupratesSI, the iron- 
based superconductors^, and other correlated electron 
system s ^ 2 ’ —. However, the previous FLEX calculations 
are restricted to the case with the spin rotational invari¬ 
ance. In order to include the SOC, we will use the tech¬ 
nique of Hugenholtz diagrams to extend the FLEX ap¬ 
proach to more general cases, where the SU(2) symmetry 
can be broken. 

Based on this method, we find that the spectral func¬ 
tion of quasiparticles is much suppressed at parts of mo¬ 
menta in the lightly doped region, suggesting the emer¬ 
gence of the pseudogap. Explicitly, the suppression oc¬ 
curs near (0,7r) point for slightly electron doping and thus 
leads to the Fermi arc near (7t/2,7t/2), while their posi¬ 
tions are reversed in the hole-doped side. These results 
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are consistent with the recent ARPES observation s^ —. 
We elaborate that the Fermi arcs and pseudogaps are 
mainly caused by the isospin fluctuation with momen¬ 
tum (7r, 7r), which overwhelms both the spin and orbital 
fluctuations. We have also studied the evolution of vari¬ 
ous collective fluctuations with doping and find that the 
isospin fluctuation dominates in the region from 30% 
hole-doping to 50% electron-doping. These results sug¬ 
gest that the scenario of J e ff = 1/2 spin-orbital Mott 
insulator is applicable to the parent and the extensively 
doped compounds of Sr 2 Ir 04 . 


II. MODEL AND METHOD 

A. Three orbital Hubbard model 


We begin with the t 2s three-orbital Hubbard model on 
the square lattice^: H = Hkin+Hsoc +Hj. The kinetic 
and SOC Hamiltonians read 


Hkin ^ 


km a. 


Hsoc = 

kmnai/3 


tisoc 


mn ' ^'a(3^ J i e7na ^kn0^ 


( 1 ) 


B. Generalized FLEX method 

In this section, we give a thorough introduction of the 
generalized FLEX method, which can naturally include 
the SOC term. The FLEX approach originates from the 
conserving approximation theory proposed by Baym and 
Kadanofl—. In this formulism, the closed-linked $ dia¬ 
grams (also known as Luttinger-Ward functional) yield 
the self-energy and the irreducible interaction vertices in 
a thermodynamic self-consistent manner, in which the 
conservation laws on the particle number, momentum, 
angular momentum and energy are respected. The FLEX 
method pioneered by Bickers and Scalapino— is the sim¬ 
plest application of the Baym-Kadanoff formulism be¬ 
yond the Hartree-Fock level, and has been widely applied 
to the single- and multi-orbital Hubbard models^,, in 
which the spin rotational invariance is respected. When 
the SU(2) spin symmetry is conserved, the scattering pro¬ 
cesses explicitly include the equal-spin particle-particle 
(PP), opposite-spin PP, equal-spin particle-hole (PH), 
and opposite-spin PH channels (see Ref—). If the SU(2) 
symmetry is broken, however, one must consider the PP 
and PH fluctuations in a comprehensive manner for there 
are mixtures of equal-spin and opposite-spin scatterings. 

We employ the Hugenholtz diagrams^ to extend this 
method to the SU(2) broken cases. Considering the on¬ 
site two-body potential 


where d' k ( 4 ma) creates (annihilates) a m-orbital 
electron with spin a and momentum k. €soc denotes 
the magnitude of SOC, and L and er are the orbital 
angular momentum operator and Pauli matrix. Ex¬ 
plicitly, the nonzero elements of L = {L X ,L V ,L Z ) for 
yz{ 1), zx( 2), and xy( 3) orbitals are L 23 = — L 32 = 
I / 31 = —L\ 3 = L z V2 = —L 21 = i. The single¬ 
particle dispersions are given by ei(fc) = —2t 3 cosk x — 
2 t 4 COS k y , e 2 (fc) = — 2 ^ 4 cos k x — 2t 5 cosk y , and e 3 (fc) = 
— 2t\[cosk x + cosky) — 4t 2 cos k x cos k y — 2t 3 (cos2fc x + 
cos 2 k y ) + /r 3 , with parameters fa, £ 3 , £ 4 , Is, p. 3 ) = 

(0.5,0.25,1.03,0.17, — 1.0)tj— Hereafter, we set £soc = 
1.03 ti without annotation, and t\ = 1 as the energy unit. 
The interaction part on the Lsite is given by 

Hi (0 fi°‘t'8pfi{U6i=j= m = n ( 1 - dap) 

ijmn ocftiiv 

+ U Si n Sjm(l ^ij') J$im$jn{^ ^ij') 

+ J &ij ^ran(1 ^ira)(l 

( 2 ) 

where U ( U ') is the intra-orbital (inter-orbital) Coulomb 
interaction, J the Hund’s coupling and J' the inter¬ 
orbital pair hopping. As usual, we set J' = J and use 
the relation U = U' + 2 J. Note that we have fabricated 
Hi to be a compact form to conveniently construct the 
Hugenholtz vertices. 


2 H (y|^l"in)cjctc„c m , (3) 

ijmn 

where the index i denote both the spin and orbital de¬ 
grees of freedom. The Hugenholtz bare vertices for the 
PP and PH channels are defined as 

r iYmn = (ij\V\mn) - ( ij\V\nm ), 

V k = r p ( 4 ) 

1 ij,mn — 1 in,jm'> 

and they are diagrammatically shown in Fig. |T] (a) and 

(b) . The Hartree-Fock (HF) and the second order d> di¬ 
agrams can be drawn easily by the bare vertices [Fig. |T| 

(c) and (d)]. Connecting the bare PP or PH vertices 
in the random-phase-approximation (RPA) series by the 
Green’s function G , we get the main body of FLEX *1) 
diagrams, which is shown in Fig. [1] (e) and (f). 

For the m-orbital system without SU(2) symmetry, the 
Green’s function and the self-energy can be expressed as 
2m x 2m matrices, satisfying the Dyson equation, 

G(k,ik n ) = [ik n 1 - h{k) - t(k,ik n )]-\ (5) 

where h(k) represents the free part Hamiltonian includ¬ 
ing SOC [see Eq. ©]• The self-energy is obtained by 
taking the derivative of $ with respect to G , i.e., pluck¬ 
ing one line of Fig. |T| (c) to (f), 
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(a) 

(b) 

(c) 





FIG. 1. (a) and (b) are bare Hugenholtz vertices for the PP 

and PH channels, (c) and (d) represent the HF term and the 
second order term, (e) and (f) are RPA forms of the closed 
linked diagrams for the PP and PH channels. The numbers 
in front of diagrams (c)-(f) are the symmetry factors. Note 
that there are no anomalous Green’s function lines because 
we only consider the normal state properties here. 


£«(*) = § £{rU^-A 0 e^° + 

q\mn 

- £(x p (g) - X p °(<?))r p ]in, jm G mn (q-k) 

+ [f h (x h (q) - X h \<i))t h ) im , jn G mn {k-q)}, 

( 6 ) 

where the abbreviation k = [k,iki] (q = [q,iqi ]) is used 
with the fermion (boson) Matsubara frequency ki = (21 + 
1)ttT (qi = 21 -kT). The first and second terms in Eq. ([6]) 
represent the HF and second order contributions, and the 
third and forth terms represent the RPA-bubbles of the 
PP and PH channels. The susceptibilities for these two 
channels are given by the (2m) 2 x (2m) 2 matrices, 

x p (q)= X p0 (g)[l + f P X p0 (q)]-\ 

X h (q)= X h0 m + f h x h0 (q)}-\ 


in which the Linhard functions are defined respectively 


by, 


Xtmn(<i) = ^£G lm (g-fc)G Jn ,(fc), 

T k (8) 

Xij,mn(q ) = £ Gim(k+q)G n j(k). 

k 

Equations J3J) to © form a closed set and thus could be 
solved self-consistently. 

The irreducible vertex of the Bethe-Salpeter (BS) 
equation for the PH channel is the differentiation of self¬ 
energy (or the second derivative of $ ) with respect to 
G. Apparently there are two kinds of contributions: (a) 
the diagrams in which the plucked two lines belong to 
the same bubble and only single RPA-form fluctuation is 
included, (b) Aslamazov-Larkin (AL) diagrams in which 
the plucked two lines belong to different bubbles and 
more than one RPA-form fluctuations are taken into ac¬ 
count. Here we omit the AL diagrams for the following 
reason s 17 ' 18 ! 28 : First, the AL contributions are demon¬ 
strated to be small (no more than 15% for the zero- 
momentum static susceptibility, for example). Second, 
the AL contributions will aggravate the agreement be¬ 
tween the FLEX results and the benchmark QMC sim¬ 
ulations. Third, the AL diagrams are necessary if we 
require keeping the conservation law rigorously, which 
is believed to be essential in some transport studies^. 
Since here we investigate the equilibrium properties of 
electrons in the normal state, the AL diagrams can be 
safely removed. Therefore we obtain the PH channel’s 
irreducible vertex 

/£ n (fc, k'-Q) = r& >mn - (t h x h (k-k')t h ] im , jn 

+ [fY(Hfc'+Q)P] in , jm , 

where momenta k, k' are for the fermions and Q is for the 
collective bosonic mode. The corresponding BS equation 
reads 

£ ll af} (k, k'-Q)G am (k , +Q)G n0 (k')^ n (k'- 1 Q) 

k ,% ,a./3 

= X h Q ^(k;Q), 

( 10 ) 

where A q and tl) h (k;Q ) represent the eigenvalue and 
eigenfunction. Particularly, if A q approaches to 1 at zero 
frequency, the system undergoes a spontaneously sym¬ 
metry breaking at momentum Q in the PH channel. 

The irreducible vertex for PP channel is the second 
derivative of $ with respect to anomalous Green’s func¬ 
tions F and , here Fij = (T T dCj) and F-- = (T t c£). 
Although there is no anomalous Green’s function lines in 
the original $ diagrams (see Fig. |T|), we can construct 
such diagrams by replacing two G with F and F t but 
keeping two arrows pointing in and other two arrows 
pointing out at each dot. Here we omit the AL diagrams 
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again and obtain the PP channel’s irreducible vertex 


spin (LS) susceptibilities are given by 


r? j)mn {k,hf\Q) = \^, m n-^ h X h {k-k')t h U, nj , ( 11 ) 

and the corresponding BS equation 


X TS (Q) a aP aX ^xt.0ma,n f ,nAQ^), 

mn ce/Hiiis 

X LS (Q) =EE vZpVpvXtp ma,nnnAQ^)> 

mn apiiu 


(13) 


- E I P ^k^'-,Q)G a Uk' + Q)G 0n {-k')r m n{k'^Q) 

k ,% ,afi 

= x Qi’!j(k\Q)- 

( 12 ) 

Unlike the PH channel, the largest eigenvalue Aq always 
associates with Q = [0,0], indicating the formation of 
Cooper pairs with opposite momenta. Eqs. HID and m 
can be used to investigate the most favorable supercon¬ 
ducting pairing gap which corresponds to tp p (k',Q) with 
the largest value of Aq. 

We have completed the introduction of the formal 
FLEX formulations, in practical application some reason¬ 
able approximations are widely used to simplify the com¬ 
putations. First, when the on-site interactions are all re¬ 
pulsive, the contributions of the PP RPA-bubbles in Fig. 
[T| (e) are relatively small and therefore could be safely left 
outd'~ 19,3 i. Second, it is more convenient to choose other 
numerical criterions, instead of Eq. m, to evaluate the 
PH channel’s instability. For example, one can choose 
the Stoner-like criterion det[i + P h X h0 {Q, 0)] < 0.003^, 
or if the biggest element of XijiQ ,0) is 50 times larger 
than Xq° (Q, 0)2*, as we employ in this paper. 

To apply the above FLEX method to our study, the 
Hugenholtz vertices in Eq. (U) are fabricated by the on¬ 
site interactions and the results presented 

in this paper are given by the parameters (U, U, J, J 1 ) = 
(5,3.5, 0.75, 0.75), no qualitative different results are ob¬ 
tained when we change the values of U and J. The nu¬ 
merical calculations are performed on 64 x 64 k meshes 
with 1024 (for T = 0.04), 2048 (for T = 0.02,0.015) and 
4096 (for T = 0.01) Matsubara frequencies. Particularly, 
we utilize the technique developed by Deisz et to ef¬ 
ficiently include the contribution of high Matsubara fre¬ 
quencies. The analytical continuation of Green’s func¬ 
tions to the real frequency is carried out by Pade ap¬ 
proximation^. The convergent solutions of the FLEX 
equations are obtained if the relative error of each ma¬ 
trix element of E is smaller than 10~ 6 . 


III. RESULTS AND DISCUSSION 

A. Collective excitations 


Let us first define corresponding susceptibilities for var¬ 
ious collective excitations relevant to the following discus¬ 
sions. The static transverse spin (TS) and longitudinal 


where the spin and orbital degrees of freedom have been 
expressed by (a, ft, fj,,v) and (m, n), respectively. The 
charge fluctuation is too small compared to other fluc¬ 
tuations and thus is not discussed here. Because of the 
introduction of SOC, the contribution of orbital fluctu¬ 
ations is no longer neglectable, so we define the static 
transverse orbital (TO) and longitudinal orbital (LO) 
susceptibilities 

X T °(Q) =£ £ L^L X mn xUc,m0n^Q^ 

afi ijmn 

(14) 

X L °(Q) LijLmnXjaiot,mPnp(Qity- 

ot.fi ijmn 


As discussed in the introduction, it has been sug¬ 
gested that the low-energy physics in Sr 2 lr 04 may be 
described as an effective one-band model with the isospin 
Jeff = 1/ 2 — i - — ■ . To check its possible application here, we 
define the isospin susceptibility. The creation operators 
for the isospin-up and isospin-down states with momen¬ 
tum k are given by 4 , ft = (4,1,4. + *4,2,+ + 4,3,t)/^ 

and 4,ft = (4,i,t - *4,2 ,t - 4,3,+)/^- where in¬ 
dexes 1, 2, 3 denote yz, zx, and xy orbitals. There¬ 
fore, the isospin operator can be constructed: = 

Efc(4+g,ft>4+q,ft)' T “( a fc,k> a fc,4) T - Then > the static 
transverse isospin (TI) and longitudinal isospin (LI) sus¬ 
ceptibilities are given by, 


X TI (Q) =§ (T T S*(q, r)S x (—q, 0 ))dr, 

X h \Q)=^ JJ (TrS z (q, T ) sz (-q,0)) d r. 


(15) 


The static TS, LS, TO and LO susceptibilities for 3% 
electron-doping at T = 0.01 are shown in Fig. [2] (a)-(d). 
As expected, the strong peaks exist around (±7r,±7r), 
which are reminiscences of the AFM order in the par¬ 
ent compound. In particular, the intensities of orbital 
susceptibilities are all stronger than those of spin sus¬ 
ceptibilities, and the TO fluctuation overwhelms the LO 
fluctuation. These features are consistent with the exper¬ 
imental results in the single-crystal neutron diffraction^ 
and nonresonant magnetic X-ray diffraction— 

Fig. [2] (e) and (f) show the TI and LI susceptibilities 
for the same 3% electron doping at T = 0.01. The peaks 
also reside at the same momenta as in the spin and orbital 
fluctuations. Quite strikingly, the intensity of the TI sus¬ 
ceptibility is almost one order of magnitude larger than 
that of all other fluctuations. To show in more detail the 
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FIG. 2. (Color online) Susceptibility for various collective 
excitations at 3% electron-doping and T = 0.01. (a), (b), 
(c) and (d) show static TS, LS, TO and LO susceptibilities, 
respectively, (e) and (f) show static TI and LI susceptibilities 
(see text). The TS fluctuation in (a) no longer exhibits the 
Ca symmetry with the introduction of SOC. 


effects of the isospin fluctuations, in Fig. [3] (a) we com¬ 
pare the maximum magnitude of x°(q) at T = 0.02, from 
40% hole-doping to 30% electron-doping, with a denoting 
TI, LI, TO, LO, TS and LS channels. The AFM order is 
determined by the numerical criterion in the PH channel, 
as introduced in Sec.II B. As shown, the TI susceptibil¬ 
ity tends to diverge as the system approaches the AFM 
state, from both the electron-doped side and the hole- 
doped side. These results suggest that the AFM order in 
the undoped and slightly doped system is mainly caused 
by the TI fluctuation. Consequently, it shows that the 
effective one-band model with the isospin J e g = l/2 LLkL£ 
can describe the magnetic properties reasonably in the 
undoped and slightly doped Sr2lr04. Another feature 
drawn from Fig. [3] (a) is that the AFM is more robust 
against the hole-doping in comparison to the electron¬ 
doping, which may be useful for the related experimental 
investigations. 

When the system is doped away from the AFM re¬ 
gion, Fig. [3] (a) shows that the TI fluctuation decreases 
rapidly with doping. However, it still has larger mag¬ 
nitude than all others in the region of 30% hole-doping 
to 50% electron-doping (the results for > 30% electron¬ 
doping are not shown here), though their differences are 
decreased with dopings. In this region, the momentum q 



FIG. 3. (Color online) (a) The largest value of various sus¬ 
ceptibilities X a ( < ?) as a function of electron doping (positive) 
and hole doping (negative) at T = 0.02, with a =TI, LI, TO, 
LO, TS or LS. (b) Fermi surfaces for 40% hole doping with 
the different colors representing the weights of the J e ff = 1/2 
doublet and J e g = 3/2 quartet, (c) The same as (b), but with 
the colors indicating the majority orbital character (magenta: 
d xz , black: d yz ). (d) The largest value of x“(q) as a function 
of the SOC magnitude £soc at T = 0.02 for 20% hole doping. 
The lines with an arrow in (b) and (c) indicate the nesting 
vector. 


is at or near the (n, 7r) point. When doping the system 
further with holes (> 30%), the spin fluctuation becomes 
dominant so that the effective one-band picture is bro¬ 
ken down. These results can be understood from the 
weights of the J e g =1/2 doublet and J e g = 3/2 quartet 
along the Fermi surface, as shown in Fig. [3] (b). One can 
see that the closed electron pocket centered around the 
T = (0,0) point is composed mainly of the J e g =1/2 
doublet, while the hole pocket around (n, n) mainly of 
the J e ff = 3/2 quartet. For 40% hole-doping, the hole 
Fermi pocket becomes large and introduces a Fermi sur¬ 
face nesting between the hole and electron pockets, as 
shown in Fig. [3] (b). Thus, the inter-pocket particle-hole 
scattering between the J e ff = 1/2 dominant band and 
Jeff = 3/2 dominant band takes action and overwhelms 
the intra-pocket scattering in the J e ff = 1/2 dominant 
band, which makes the isospin fluctuation no longer the 
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FIG. 4. (Color online) Spectral functions A(k, 0) for 10% 
electron-doping (a) and 5% electron-doping (b) at T = 0.015. 
Spectral functions A(k, 0) at T = 0.02 (c) and T = 0.01 (d) 
for 3% electron-doping. The pink (black), green (red), gray 
(blue) lines in (e) are u> dependence of the spectral functions 
for T = 0.02 (T = 0.01), at three momenta indicated by 
the arrow in (c) and (d). An enlarged version for T = 0.01 
is shown in (f). The white lines in (c) and (d) denote the 
boundary of the magnetic Brillouin zone. 


leading collective excitation. On the other hand, from 
the distribution of orbital characters shown in Fig. [3] (c), 
we can find that the main inter-pocket scattering occurs 
in the same orbital, which makes the orbital fluctuations 
also suppressed. Therefore, the spin fluctuation is the 
dominant collective excitation for the heavily hole-doped 
system. However, the Fermi level moves away from the 
Jeff = 3/2 dominant band with the electron doping, con¬ 
sequently it is only the J e g =1/2 dominant band that 
crosses the Fermi level. Therefore, the one-band picture 
is always robust against the electron doping. Further¬ 
more, we have also examined the range of validity of the 
effective one-band picture with respect to the strength of 
SOC. As shown in Fig. 0(d), the TI susceptibility prevails 
all others for 20% hole-doping if £,soc > 0.8, indicating 
that this picture survives in an extended range. 


B. Weak pseudogap behavior and Fermi arcs 

The pseudogap behavior can be detected from the 
single particle spectral function which is defined as 
A(k,ui) = — ImG nn (fe,u;)/7r, where n denotes both 
spin and orbital indices. In Fig. 0 (a) and (b), we 
show the contour plot of the zero-energy spectrum at 
T = 0.015 for 10% and 5% electron-doping, respectively. 
Since the system has been away from the AFM order at 
these dopings [see Fig. 0 (a)], an intact Fermi surface 
is expected in the conventional normal state. Indeed, 
for 10% electron-doping, a closed diamond-shaped Fermi 
surface is observed, as shown in Fig. 0(a). However, an 
obvious reduction in the spectral intensity occurs around 
(0,7r) and its symmetric points for 5% electron-doping. 
This reduction is strongly temperature-dependent, be¬ 
cause it appears only below a certain temperature, as 
shown in Fig. 0 (c) and (d) at 3% electron-doping for 
T = 0.02 and T = 0.01, respectively. Owing to the 
suppression of the spectra near (0,±7r) and (±7r,0), the 
Fermi surface around these momenta is destructed and 
the four residual separated segments form the so-called 
Fermi arcs [see Fig. 0 (b) and (d)]. These results are 
consistent with the recent ARPES experiment^. In or¬ 
der to look in more detail the suppression of the spec¬ 
tra, we show the energy distribution curves (EDC) of the 
spectral functions for 3% electron-doping at three mo¬ 
mentum points near (7r, 0) in Fig. 0 (e). One can see 
that the suppression occurs basically around the Fermi 
energy when the temperature is decreased from T = 0.02 
to T = 0.01. Moreover, different from a well defined 
quasiparticle peak at T = 0.02, the spectral functions 
at T = 0.01 for q = (307r/32,0) and (3l7r/32,0) show a 
weak dip around the Fermi energy, as shown in Fig. 0 
(f). It suggests that a weak pseudogap does exist around 
(7 r, 0) and its symmetric momenta for slightly electron¬ 
doping. 

Since thers is no long-range order in this doping range, 
the pseudogap is most likely to result from the scattering 
of quasiparticles by the collective excitations. In this 
framework, the quasiparticles around (7c, 0) and (0,7r) 
are strongly coupled by the collective excitations with 
the transferred momentum (7r, 7r) in the scattering pro¬ 
cess. In accord with this analysis, the results presented 
in Fig. 0 demonstrate that all collective excitations, in¬ 
cluding the spin, orbital and isospin fluctuations, exhibit 
the characteristic momentum (7r,7r). In particular, the 
transverse isospin (TI) fluctuation overwhelms the spin 
and orbital fluctuations. We have also carried out a care¬ 
ful examination, and find that all the important elements 
of x h (q, 0) are included in the TI fluctuation. Especially, 
for q = (7r, 7r) the elements included in the TI susceptibil¬ 
ity are calculated to amount to 91.3% weight of the total 
PH susceptibilities for 3% electron-doping at T = 0.01. 
Therefore, it suggests that the scatterings of quasiprticles 
by the (7r, 7r) TI fluctuation lead to a weak pesudogap 
partially opened around (0,±7r) and (±7r, 0). 

We have also investigated the possible pseudogap be- 
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FIG. 5. (Color online) (a) and (b) are spectral functions 
A(k, 0) at T = 0.04 and T = 0.01 for 17% hole-doping. The 
white lines denote the boundary of the magnetic Brillouin 
zone. 

havior and Fermi arcs in the lightly hole-doped region. 
The typical results of the spectral function A(k,0) for 
17% hole-doping are shown in Fig. [5] At T = 0.04, 
there is a nearly circular Fermi surface around the T 
point. By decreasing the temperature to T = 0.01, the 
spectral weights are suppressed on some Fermi momenta 
and consequently it leads to the formation of the Fermi 
arcs. However, the suppressions now appear around 
(7r/2,7r/2), which is contrary to the case of the electron¬ 
doping shown in Fig.[]](b) and (d). This suppression also 
results from the (7r, 7r) TI fluctuation. In this case, the 
Fermi surface shrinks in comparison with that of electron¬ 
doping, so that the “hot spots” (the crossing points of 
the Fermi surface with the boundary of the AFM Bril¬ 
louin zone) at which the quasiprticles are strongly scat¬ 
tered move from the (±7r,0) and (0, ±7r) points to the 


(±7t/2, ±7t/2) points (see Fig. [4] and Fig. [5j). 


IV. SUMMARY 

In summary, we have extended the FLEX approach 
by Hugenholtz diagrams to include the case where SU(2) 
symmetry is broken. Using this approach, we investigate 
various collective fluctuations and the spectral function 
of quasiparticles. It is found that the isospin fluctuations 
derived from the J e ff = 1/2 spin-orbit Mott insulator 
dominate over the spin, orbital and charge fluctuations 
in the extended doping region, suggesting the validity 
of the isospin J e g = 1/2 picture in an extended doping 
regime. Also the isospin fluctuation leads to the emer¬ 
gence of the pseudogap and Fermi arcs in the slightly 
doped system, which is consistent with the ARPES ex¬ 
periments for slightly doped S^IrO^dil. 
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